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FOREWORD 


This  report  summarizes  the  research  conducted  by  the  Georgia  Institute  of 
Technology  (GIT)  Center  of  Excellence  in  Rotary  Wing  Technology  (CERWAT)  during 
the  second  phase  of  the  Army's  rotorcraft  Centers  of  Excellence  (COE's)  program.  The 
major  accomplishments  achieved  during  the  first  two  phases  of  CERWAT  are  illustrated  in 
Figure  1  along  with  some  items  addressing  the  scope  of  Phase  ni,  the  Program 
Sustainment  Phase,  which  began  on  1  November  1992  and  runs  for  three  years.  You  will 
note  that  for  this  third  phase,  the  name  CERWAT  has  been  changed  to  CERT,  Center  of 
Excellence  in  Rotorcraft  Technology,  to  reflect  the  Army's  name  for  Phase  ni  of  the 
rotorcraft  COE  program. 

As  can  be  seen  in  Figure  1,  there  were  five  key  objectives  identified  for  Phase  n  -  ^  _ . 

the  Program  Maturation  Phase.  The  first  objective  was  to  add  the  flight  mechanics  and 
controls  discipline  to  the  research  program.  This  addition  provided  the  missing  critical 
discipline  for  advancing  rotorcraft  technology  and  served  as  a  catalyst  for  integrating 
controls  technology  with  the  other  rotorcraft  disciplines  -  aerodynamics,  aeroelasticity,  and 
structures  and  materials.  The  second  objective  was  to  en^hasize  interdisciplinary  research 
among  the  four  critical  rotorcraft  disciplines.  This  has  been  accon:Q}lished  to  a  large  extent 
as  discussed  in  the  write-ups  on  the  individual  research  tasks.  The  third  objective  was  to 
stress  rotorcraft  related  research  development.  The  negotiated  funding  profile  for 
CERWAT  is  illustrated  in  Figure  2.  As  can  be  seen,  the  funding  was  greatly  reduced  in 
FY92  and  there  was  no  assurance  of  a  Phase  m  program.  As  a  result  of  active  research 
program  development  by  all  of  the  CERWAT  researchers,  rotorcraft  related  research  has 
grown  substantially  during  the  Phase  n  program  as  illustrated  by  Figure  m,  although  that 
growth  has  leveled  off  and  rotorcraft  related  research  is  expected  to  diminish  in  FY93  due 
to  Department  of  Defense  funding  cuts.  It  would  be  difficult  to  sustain  the  GIT  rotorcraft 
COE  without  an  Army  Phase  III  Program.  The  fourth  objective  for  Phase  n  was  to 
emph^ize  Army  officer  graduate  student  recruitment.  During  CERWAT  Phase  n,  14 
Army  officers  received  M.S.  or  Ph.D.  degrees.  This  was  accomplished  through  talks 
given  at  the  Army  Aviation  Officers  Advanced  Course  at  Fort  Rucker,  AL,  and  mostly 
through  word  of  mouth  firom  Army  officers  who  graduated  and  knew  the  quality  of  the 
CERWAT  program.  The  final  objective  of  CERWAT  Phase  II  was  to  emphasize 
technology  transfer  to  industry  and  government  the  research  results  that  were  ready  for 
application  in  the  next  phase  of  the  research  and  development  cycle.  The  Georgia  Tech 
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CERWAT  has  been  very  active  in  this  area  and  specific  discipline  results  are  addressed  in 
the  task  write-ups.  Some  overall  examples  are  provided  in  Hgure  4. 

In  summary,  it  can  be  seen  that  the  Georgia  Tech  CERWAT  accomplished  the 
objectives  it  set  for  itself  during  the  Phase  n  program.  We  feel  that  we  have  a  mature 
program  and  look  forward  to  executing  the  Phase  HI  -  Program  Sustainment  Phase. 
Objectives  we  have  established  for  Phase  n  are  as  follows; 

•  Work  with  industry  and  government  in  identifying  needs  and  areas 
where  CERT  can  contribute 

•  Build  on  individual  disciplinary,  multi-disciplinary,  and  interdisciplinary 
successes 

•  Explore  opportunities  available  as  a  result  of  new  rotorcraft  and  related 
initiatives 

•  Play  more  of  a  leadership  role  in  advancing  rotorcraft  technology. 


osoimiiAinzcN 


CERWAT  PROGRAM  PHASES 

PHASE  I  •  PROGRAM  PSTIATION  AND  DEVELOPMENT  (1Ma-19t7)  CERWAT 

•  CURRICULUM  DEVELOPMENT  (3  COURSES -DISCOURSES) 

•  FAOUTY  UPGRADE  AND  DEVELOPMENT  (2  FACILinES  UPGRADED, 
a  NEW  FAaunES  DEVELOPED) 

•  FACULTY  RECRUITMENT  (i  FACULTY  *14  FACULTY) 

•  CERWAT  FELL0W8MP  PROGRAM  DEVELOPMENT  (IS  FELLOWS) 

•  RESEARCH  PROGRAM  DEVELOPMENT  (3  DBOPUNES:  AERODYNAMICS, 
AEROELASnCfTY,  AND  STRUCTURES  «  MATERIALS) 

PHASE  S  -  PROGRAM  MATURATION  (ISSS  •  1SBS)  CERWAT 

•  FLIGHT  MECHANICS  A  CONTROL  DSCIPUItt  ADDED 

•  SfTERDSCWLmARY  RESEARCH  EMPHASOEO 

•  ROTORCRAFT  RELATED  RESEARCH  DEVELOPMENT 

•  ARMY  OFFICER  GRADUATE  STUDENT  RECRUITMENT 

•  TECHNOLOGY  TRANSFER  EMPHASIZED 
PHASE  M  .  PROGRAM  SUSTASIMENT  (1893-1SSB)  CERT 

•  ARMY  FUNDING  CUT  M  HALF 

PIDUSTRY  COST  SHAIUNG  INQUIRED  AND  UNIVERSITY  COST  SHARING 
MCREASED 

•  A  FLIGHT  SIMULATION  TASK  ADDED 


Hgure  1. 
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QIOgNMIATECIKi  . . .  . 

CERWAT  miTlATIVES  TO  EXPAND 
ROTORCRAFT  RESEARCH  ACTIVITIES 
FOR 

ADVANCING  ROTARY  WING 
TECHNOLOGY 

•  Sbmilatlon  and  MoMIng 

•  Flight  Simulation  (FLIGHT  9M)  Laboratory 

•  2GCHAS  Enhancamant  Program 

•  PM-TRADE/DARPAADST  Program 

•  Artificial  Intalliganca  and  Robotica 

•  RPA  Miaaion  EffacHvanaaa 

•  Aarlal  Robotica  Competition 

•  Low  Obaarvablea 

•  CFO  Modeling 

•  SafMy  Enhancamenta 

•  USAF  PAVE  LOW/PAVE  HAWK  Structural  bdagrity  Program  (SIP) 

•  TrI'Sarvlea  Flight  Data  Recorder  Program 

•  DhreraMcatlon 

•  Tlltrotor  and  HSRC 

•  FAA  and  Commercial  Proqrama. 


Hgure4. 
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Khan,  Azhar,  Rotor  Performance  and  Helicopter  Mission  Optimization  with  a  Quasi-Static 
Inflow  Theory,  Ph.D.  lliesis,  Georgia  Institute  of  Technology,  December  1992. 

Li,  J.,  Inter  .aminar  Fracture  Analysis  for  Lamiruited  Composites  under  Combined 
Loading,  PhJ3.  Thesis,  Georgia  Institute  of  Technology,  September  1992. 


M.S.  Special  Problems: 

Foley,  Susan  M.,  "The  Correlation  Between  Canard  Position,  Voncx  Trajectory  and 
Pressure  Distribution  Over  a  Wake  /  Wing  /  Canard  Model".  M.S.  Special  Problem  Report, 
School  of  Aerospace  En^eering,  Georgia  Institute  of  Technology,  September  1992. 

Petroski,  Matthew,  "Swept  Tip  Blade  Flow  Analysis,  Vortex  Trajectory  Mapping,  and 
Laser  Light  Intensity  Study".  Special  Problem  Report,  School  of  Aerospace  Engineering, 
Georgia  Institute  of  Technology,  September  1992. 


Abstracts  Submitted/Accepted 

Affes,  H.,  Conlisk,  T.,  Kim,  J.M.,  and  Komerath,  N.M.,  "Computation  and  Experiments 
on  the  Pre-Collision  Phase  of  a  Vortex  Interacting  with  a  Cylindtf .  Accepted  by  the  AlAA 
24th  Huid  Dynamics  Conference,  Orlando,  FI. 

Armanios,  E.  A.,  Hooke,  D.,  Kamat,  M.,  Palmer,  D.,  and  Li,  J.  "Design  and  Testing  of 
Composite  Laminates  for  Optimum  ^tension-Twist  Coupling,"  Submitt^  for  publication 
in  an  ASIM  STP  on  Con^XKite  Materials:  Testing  and  Design. 

Byms,  E.V.  Jr.,  Calise,  A.J.,  'Turtlier  Development  of  an  Approximate  Loop  Transfer 
Recovery  Methodology  for  Fixed  Order  Dynamic  Compensator  Design,"  Optirnal  Control 
^plications  and  Methods. 

Byms,  E.V.  Jr.,  Calise,  A.J.,  "Approximate  Recovery  of  Hoe  Loop  Shapes  Using  Fixed 
Order  Dynamic  Cori^nsation,"  AIAA  /.  ofGuidcuice,  Control,  and  Dynamics. 

Fawcett,  P.A.,  Komerath,  N.M.,  "Spatial  Cross-Conelation  Velocimetry:  Application  to 
Unsteady  Flows".  Submitted  to  Journal  ofAircrcfl. 

Funk,  R.,  Fawcett,  PA.,  and  Komerath,  N.M.,  "Instantaneous  Velocity  Fields  in  a  Rotor 
Wake  By  Spatial  Correlation  Velocimetry".  Accepted  by  the  AIAA  24th  Fluid  Dynamics 
Conference,  Orlando,  R. 

Kim,  J-M.,  and  Komerath,  N.M.,  "Overview  oK  the  interaction  between  a  rotor  vortex 
system  and  an  airframe".  Accepted  by  the  AIAA  24th  Fluid  Dynamics  Conference, 
Orlando,  FI. 

Komerath,  N.M.,  Fawcett,  P.A.,  "Spatial  Cross-Correlation  Velocimetry:  Theoretical 
Basis  and  Validation".  Submitted  to  Journal  cfAircrcfl. 

Kim,  J-M.,  Komerath,  N.M.,  and  Liou,  S-G.,  "Vorticity  Concentration  at  the  Edge  of  the 
Inboard  Vortex  Sheet".  Submitted  to  AHS  Joumcd.. 
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Kim,  J-M.,  Komerath,  N.M,  and  Liou,  S-G.,  "Vorticity  Concentration  at  the  Edge  of  the 
Inboard  Vortex  Sheet".  Accepted  for  presentation  at  the  49th  AHS  Forum,  Sl  Louis,  MO, 
May  1993. 


Lai,  M.K.,  Liou,  S-G.,  Pierce,  G.A.,  and  Komerath,  N.M.,  "Measurements  Around  a 
Rotor  Blade  Excited  in  Pitch.  Part  2:  Blade  Surface  Pressure".  Submitted  to /owvw/o/fAe 
American  Helicopter  Society. 

Lai,  M.K.,  Liou,  S.G.,  Pierce,  G.A.,  Komerath,  N.M.,  "Measurements  of  Unsteady 
Pressure  on  a  Pitching  Rotor  Blade".  Accepted  for  inclusion  in  the  Proceedings  of  the  49th 
AHS  Forum. 


Lai,  M.K.,  Liou,  S-G.,  Komerath,  N.M.,  and  Pierce,  A-G.,  "Correlations  of  Pressure  and 
Velocity  Data  with  predictions  for  Unsteady  Inflow  to  a  Vibrating  Rotor  Blade".  Accepted 
for  presentation  at  Ae  AIAA  24tii  Fluid  Dynamics  Conference,  Orlando,  FL 

Leitner,  J.  and  Prasad,  J. V.R.,  "Helicopter  Robust  Nonlinear  Controller  Synthesis  using 
Approximate  Linearization,"  Paper  accepted  for  presentation  at  the  49th  Annual  Forum  of 
the  American  Helicopter  Society,  St.  Louis,  1993. 

Liou,  S-G.,  Komerath,  N.M.,  Lai,  M.K.,  "Measurements  Around  a  Rotor  Blade  Excited 
in  Pitch.  Part  1:  Dynamic  Inflow".  Submitted  to  Journal  of  the  American  Helicopter 
Society. 

Pamas,  L.,  Armanios,  E.  A.,  Sriram,  P.  and  Rehfield,  L.,  "Analysis  and  Testing  of 
Postbuckling  and  Criroling  in  Advanced  Composite  Stiffeners,"  Accepted  for  publication. 
Journal  of  Aerospace  Engineering,  ASCE,  1993. 

Petroskd,  M.,  Liou,  S-G.,  Komerath,  N.M,,  "Ccraparison  of  the  Flow  over  a  Blade  Tip  at 
High  Pitch  Angles  under  Rotating  and  Fixed-Wing  Conditions".  Accepted  for  presentation 
at  tile  AIAA  24th  Fluid  Dynamics  Conference,  0:^do,  FL 
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8.  SCIENTIFIC  PERSONNEL  SUPPORTED  BY  THIS  PROTECT  AND 
DEGREES  AWARDED  DURING  THIS  REPORTING  PERIOD; 


Co*Principal  Investigators:  R.  B.  Gray  and  D.  P.  Schrage 

Faculty:  R.B.  Gray,  D.P.  Schrage,  E.A.  Armanios,  A.J.  Calise, 

J.I.  Craig,  S.V.  Hanagud,  D.H.  Hodges,  N.M.  Komerath, 
D.A.  Peters,  G.A.  Pierce,  J.V.R.  Prasad,  L.N.  Sankar, 
J.C.  Wu 

Research  Engineers:  C.M,  McKeithan,  C.M.  Wang 


Post  Doctoral  Fellows:  S.G.  Liou  ,  and  P.  Srirani 


Graduate  Research  Assistants: 

Ph.D.:  Y.  Chung,  P.  Fawcett,  R.  Funk,  D.  Gillam,  L.  Gunomadi, 

N.  Hariharan  J.  Higman,  D.  Hixon,  G.  Kim,  J.  Kim,  M.  Lai 
J.  Li,  J.  Lcimer,  M.  Mittal,  D.  Palmer,  W.  Park,  M.T.  Patterson, 
D.  Reece,  S.  Savanur,  X.  Shang,  W.  Stumpf,  F.L.  Tsung, 
N.  Weston,  C.  Won,  X.  Xiao,  J.  Zhang. 

M.S,:  J.  St.-Angelo,  S.  Foley,  M.K.  Lai,  G.  Nichols,  M.  Petroski, 

and  B.  Zhao 


CONTRIBUTED  TO  PROJECT,  BUT  WERE  NOT  SUPPORTED 
Faculty: 

Prof.  V.L.  Berdichevsky 
Post  Docs: 

A.  Atilgan  V.  Sutyrin, 

Ph.D.  Students: 

C.E.S.  Cesnik,  M.  Cibulka,  Major  W.  Lewis,  and  T.  Smith, 

M.  S.  Students: 

T.  Bowles,  D.  Boyd,  P.W.  Cheang,  Captain  S.  Ingalls, 

B. D.  Palmer,  ILt  D.  Troutman,  and  Captain  D.  Wheelock 
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DEGREES  AWARDED  (THIS  REPORTING  PERIOD) 


NAME 

DEGREE 

-  DATE 

LAST  KNOWN  AFFILUTION 

D.  de  Andrade 

Ph.D. 

SEP  92 

Post  Doc  >  Georgia  Tech 

P.A.  Fawcett 

Ph.D. 

SEP  92 

S.M.  Foley 

M.S. 

SEP  92 

J.  Leitner 

M.S. 

SEP  92 

Ph.D.  Program,  Ga.  Tech 

W.D.  Lewis 

Ph.D. 

SEP  92 

AATD,  Fl  Eusds,  VA 

D.S.  Reece 

M.S. 

SEP  92 

Ph.D.  Program,  Ga.  Tech 

W.M,  Stumpf 

Ph.D. 

SEP  92 

Post  Doc  -  Georgia  Tech 

Y.  Chung 

M.S. 

DEC  92 

Ph.D.  Program,  Ga.  Tech 

M.  Petroski 

M.S. 

DEC  92 

PhJD.  Program,  Ga.  Tech 

L.N.  Gummadi 

Ph.D. 

DEC  92 

Post  Doc.  Georgia  Tech 

J.Li 

Ph.D. 

DEC  92 
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RESEARCH  TASKS 


L  Aerodynamics 


Task  1.  Studies  in  Three-Dimensional  Viscous  Aerodynamics 
J.C.  Wu  and  CM  Wang 
Problem  Studied 

A  numerical  approach  based  on  vorticity-velocity  fonnulation  and  an  integral 
representation  of  the  velocity  vector  has  been  develop^  in  the  past  by  the  present 
researchers  for  the  study  of  complicated  flow  problems.  Many  2-D  flows  including 
dynamic  stall  of  airfoils  can  be  accuratel^ync^iicted  by  the  approach.  The  present  phase 
involves  two  distina  efforts.  The  first  ef^m  focuses  (» the  extension  of  the  numerical 
approach  to  3-D  incoopessible  flows.  The  second  effori  is  die  development  of  the  mediod 
for  compressible  flows. 

Progress  During  the  Last  Reporting  Period 

The  extension  of  the  present  method  is  hot  staightforward.  The  major  difficulties 
include:  1.  accurate  determination  of  the  vordcity  on  solid  surfaces:  2.  satmacticm  of  the 
divergence-free  condition  for  the  vordcity  field  (where  the  vorticity  vector  has  duee 
disdncdve  components);  and  3.  accurate  detemdnadtm  of  the  surface  pressure.  These 
issues  have  been  successfully  resolved  in  the  present  phase.  Flows  around  flat-plate  and 
rectangular  wings  of  NACA  (X)12  secdon  have  been  computed  by  the  attended  method. 
Soludon  of  the  flat-plate  wing  has  been  presented  in  an  AIAA  p^er  (91-3262)  and  will  not 
be  included  here.  Soludons  of  the  wug  of  NACA  0012  secdon  with  rounded  dp  are 
briefly  discussed  as  follows. 

Flows  of  three  angles  of  attack  (S^  15**  aiui  25*0  have  been  cooqruted  on  a  rather 
coarse  0-0  grid  of  41  X  12  X  31,  whm  the  three  numbers  denote  the  grid  number  in 
choidwise,  spanwise,  and  normal  direcdtm  respecdvely.  For  die  angle  of  attack  S”  case, 
flow  is  tot^y  attach^  and  the  dp  vortex  is  ramer  weak  udiere  the  vortex  is  attached  and 
wrapped  around  the  wing  dp  before  it  leaves  the  tndling  edge.  The  dp  vortex  becomes 
stronger  for  the  15**  case,  but  the  vortex  is  sdll  attached  fully  mi  the  wing  dp.  However, 
large  choidwise  flow  separadon  is  observed  on  the  uroer  surface  and  tte  reversed  flow 
region  is  dindiiishing  toward  the  The  particle  trace  plot  near  the  upper  surface  indicates 
an  eye-type  pattern  u  formed  which  has  uso  been  observed  in  wiixl  tunnel  testing.  Fordie 
25°  case,  the  strong  tip-vortex  lifts  off  from  the  tipnear  the  60%  chord  and  sectitmal  lift 
shows  a  hunqi  near  tte  tip  caused  by  the  lift-off.  Flow  sqiantes  from  the  wing  leading- 
edge  for  most  of  the  span. 

The  compressibility  effea  is  addressed  in  die  second  effort  of  die  present  phase. 
This  effect  is  accounted  for  by  the  use  of  decomposition  of  the  velocity  vector  in  a 
Helmholtz  manna*.  The  rotational  part  of  die  vefoci^  can  be  solved  by  the  present 
incompressible  vorticity-velocity  appre^h  uimeas  the  inotadonal  part  is  omy  of  limited 
size  (comparable  to  die  vortical  n^on),  compared  to  the  solution  domam  of  a  full 
compressible  flow  solver.  The  ^bridizadon  of  die  incooTressible  Navier-Stdees  solver 
and  tlw  modified  Euler  Sdver  Ofifen  advantages  in  computational  effkiency  and  accuracy. 
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Numerical  results  of  flows  around  circular  cylinders  and  airfoils  show  excellent  agreement 
with  experimental  results  over  a  rather  broad  Mach  number  range.  However,  no  strong 
shocks  can  be  treated  at  die  present  time. 

The  detailed  mathematical  and  numerical  foimulatitms  of  die  present  efforts  as  well 
as  the  numerical  results  and  discussions  will  be  included  and  elaborated  upon  in  the 
following  upcoming  Ph.D.  theses: 

1 .  Kim,  G.,  A  Vorticity-Velocity  Approach  for  Three-Dimensional 
Unsteady  Viscous  Flows  over  Wings,  Ph.D.  Thesis,  Georgia  Institute  of 
Technology,  in  preparation. 

2.  ^g,  Q.,  Confutation  of  Flows  Around  Oscillating  and  Fast  Pitched 
Flat  Plates,  Phi).  Thesis,  Geoigia  ^titute  of  Technology,  in  preparatitHi. 

3.  Patterson,  M.T.,  An  Integro-Differential  Solution  Methodology  for 
Compressible  Rotational  Flows,  PhD.  Thesis,  Georgia  Institute  of  Technology,  in 
preparation. 
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Task  2.  Aerodynamic  Interactions  * 

P.L:  Dr.  N.M.  Komerath 

Research  Team:  Robert  Funk,  Jai-Moo  Kim,  Shiuh-Guang  Liou 


Obiectivcs 

This  task  seeks  to  advance  ci^abilities  and  knowledge  in  the  area  of  aerodynamic  inteiac- 
titms  in  rotorcraft  flows. 

Progress  duriny  the  final  period!  7/92.1/93 

The  effort  during  this  period  was  concentrated  on  vortex-surface  interactions  and 
vortex-shear  layer  interactions.  Jai-Moo  Kim's  dissertation  is  to  be  defended  in  mid- 
February.  Papers  have  been  accepted  for  ^sentation  at  the  AIAA  Fluid  Mechanics 
conference  on  vortex-surface  interactions  (joint  woric  widi  Prof.  Conlisk's  group  at  OSU), 
spatial  conelation  velocimetry  in  the  rotor  wake,  and  summarizing  the  behavior  of  the  rotor 
wake  around  the  airframe.  The  results  obtained  in  this  period  are  summarized  along  with 
the  final  contract  summary  below. 

Summary  of  achievements  duriny  the  5-vcar  contract  period 

1.  Rotor  Wake  /  Airframe  Interaction 

A  good  definition  has  been  obtained  of  the  interaction  between  the  vortex- 
dominated  wake  of  a  rotm- and  the  flow  over  an  airfiame.  For  the  case  where  there  is  no 
la^e-scale  separation  of  die  mean  flow,  a  final  piece  of  the  puzzle  came  recendy  whoi  J.M. 
Kim  showed  the  impc^ant  role  of  the  cote  axial  velocity  in  explaining  die  asymmetric 
vortex  behavior,  and  discovered  the  remains  of  the  tip  vortex  below  the  airframe.  Now  we 
can  provide  guidance  to  those  computing  die  wake/airfiame  interaction  cm  what  to  do  with 
Ae  vortex  near  die  airframe. 

The  presence  of  separated  mean  flow  was  found  to  be  surpti^gly  sinqile  to  model 
(Kim,  thesis  *93).  At  the  low  advance  ratios  where  wake/fuselage  interactitm  is  smmg,  it 
was  found  that  the  vorticity  in  the  sqiarated  shear  layer  was  relatively  weak,  so  that  the  tip 
vortex  totally  dominated  the  flow  tehavior.  It  appears  that  the  separated  regimi  can  be 
reasonably  modeled  by  a  weak  recirculation  zone  bounded  by  a  vortex  sheet  The  tip  vortex 
suffers  little  modification  in  passing  through  this  zone,  and  the  vortex/airframe  interaction 
proceeds  in  a  similar  fashion  to  the  case  without  mean  flow  sqiaration.  The  separation 
zone,  on  the  other  hand,  is  periodically  destroyed  and  re-established.  We  anticipate  that  the 
dominant  features  the  surface jmssuie  can  be  c^tured  by  sinqile  models  without  reson 
to  Navier-Stokes  formulations.  Of  course,  there  are  several  secmidary  features  that  need 
further  examination,  especially  widi  regard  tt>  the  possibility  of  flow  control  exploitation. 

2.  Vortex-surface  interaction 

The  intmction  of  a  vortex  with  a  curved  surface  lemaiiis  as  die  dominant  problem 
needing  detailed  experimentation  and  arudysis.  Egress  was  made  in  achieving  excellent 
correlation  between  calculation  (OSU)  and  eiqieriinents  up  to  the  stage  where  die  vortnc 
core  started  interacting  with  the  surface  boundary  layer.  In  die  last  tqiorting  period,  a 
concerted  effort  was  made,  with  nof.  Conlisk  present  to  provide  guidance,  to  vwialize  ^ 
initiation  of  boundary  layer  separation  upstream  of  the  vortex.  After  2  days  of  intensive 
efforts,  several  frames  of  vi^  were  obtained  usin^  an  intensified  camera  showing  clear 
separation  patterns,  ct^itured  widi  a  tenyoral  resolution  of  2S  nanoseconds. 
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3.  Rotor  Wake  GetHoetry 

Kim's  finding  of  the  concentrated  roll-up  of  the  inboard  vortex  sheet  explained 
some  of  the  "secondary"  vortex  features  previously  observed  by  Liou  and  Brand.  Kim's 
subsequent  experiments  clearly  showed  that  die  roll-up  was  gready  accelerated  by 
interaction  of  die  sheet  with  the  tip  vortex  from  the  following  blade.  In  fact,  the  subsequent 
vortex  trajectory,  and  thus  the  "skewing"  and  "contraction"  of  the  near  wake,  are  bound  to 
be  influenced  by  this  counter-rotating  vortex  to  a  higher  degree  than  by  mutual  interaction 
between  tip  vortices.  This  was  an  unexpected  finding,  which  originate  in  a  disagreement 
between  a  Ph.D.  candidate  and  his  advisor,  followed  by  the  candidate  undertaking  his  own 
experiments  to  prove  his  point  beyond  doubt 

4.  Spatial  Gonelatitxi  Velocimetry  ' 

The  high-speed  planar  velocimetry  effort,  undertaken  as  a  high-risk  exploration  in 
1988,  has  succeeded,  and  is  producing  a  unique  capability.  An  unexpected  bonus  is  the 
feasibility  of  application  to  large-area  measurement  without  lasers.  Also,  since  the 
technique  works  without  phase-averaging,  we  see  the  way  open  to  undertake  much  more 
challenging  explorations  of  turbulent  phenomena.  During  the  final  period,  Philip  Fawcett 
succeeded  m  capturing  a  larjge  number  of  instantaneous  velocity  fields  in  the  wake  of  an 
isolated  rotor  in  forward  flight  A  paper  based  on  these  results  has  been  accepted  for 
presentation  at  the  Fluids  &  Plasma  Conference  in  July  '93.  A  U.S.  Patent  application 
based  on  this  technique  has  been  approved,  and  is  in  the  ^al  stages  of  revision  of  figures. 

5.  Multi-Diagnostic  Capability 

The  phase-resolved  point  measurement  techniques,  the  desynchronized 
visualization  technique,  and  SCV  have  come  together  in  a  multi-faceted  measurement 
capability,  which  opens  the  way  for  experinrental  research  on  complex  configurations 
during  maneuvers.  The  first  demonstration  of  the  power  of  tius  approach  came  in  our  study 
of  wake/wing  interaction  (Susan  Foley,  M.S.  Special  Problem  and  AIAA  92-4(X)8),  which 
revealed  unexpected  phenomena  relat^  to  flow  separation  and  reattachment  WMle  the 
individual  pieces  may  appear  mundane,  the  order-of-magnitude  increase  in  tire  productivity 
of  advanc^  diagnostic  techniques  has  profound  implications  for  the  feasibility  of  studying 
many  interaction  problems. 
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Task  3.  Blade  Tip  Aerodynamics 


Probiem  Studied 

The  objective  of  this  research  is  to  understand  the  behavior  of  the  flow  over  rotor 
tip  shapes  for  use  at  high  pitch  angles,  and  to  develop  methods  of  improving  retreating- 
blade  performance  in  forward  flight 


During  the  past  five  years  of  research  effort  much  progress  was  made  under  this  task.  The 
work  performed  under  this  project  has  been  documented  in  open  literature.  Several  joint 
studies  and  publications  with  researchers  in  Ae  U.  S.  Hfelicopter  industries  were  also 
carried  out  during  this  period.  A  list  of  all  publications  under  GERWAT  support  for  the 
period  October  1987-  October  1992  is  enclosed. 

Some  of  the  major  accomplishments  during  this  period  are  as  follows: 

1)  The  3-D  Navier-Stokes  solver  developed  under  the  origiiud  CERWAT  program  has 
bMn  continually  updated  to  in^nrove  its  speed,  and  accur^.  We  have  replaced  die  second 
order  accurate  miite  difference  formulas  in  this  solver  widi  fourth  order  accurate  operator 
compact  forms.  The  original  numerical  viscosity  terms  in  this  solver  have  been  repUoBd 
with  a  scheme  based  on  the  Roe  upwind  scheme.  The  Baldwin-Lomax  turbulence  model 
has  been  replaced  with  a  k-e  turbulence  nxxleL  For  hover  and  quasi-steady  forward  flight 
studies,  the  CPU  time  needed  has  been  reduced  to  less  than  1  hew  on  the  (Say  Y/MP  class 
of  cooqtuters  (for  a  120,(X)0  point  grid)  tiirough  local  time  stepping  and  r^lacement  (tf  tiie 
block  tri-diagonal  matrix  inversions  within  dte  flow  solver  wim  scalar  tri-diagonal  nutrix 
inversions. 

2)  The  3-D  Navier-Stokes  solver  was  applied  to  advanced  rotor  blade  tip  shapes.  A 
Georgia  Tech  ^  shape,  and  a  BERP-like  tip  were  studied.  Where  possible,  conqtarisons 
of  the  velocity  field  were  made  with  the  measurenaents  carried  out  by  Ptof.  Krxnerath  We 
also  studied  me  effects  of  blade  rotation  on  the  flow  over  the  blade  tip,  by  coixq)uting  the 
flow  over  the  BERP-like  planfotm  in  the  fixed  wing  mode  and  rotar^r  wing  mode.  It  was 
found  that  the  cenoifiigal  flow  effects  in  the  rotary  wing  noode  had  a  significant  beneficial 
effect  on  the  rotary  patonoance  particularly  tX  hi^  coU^ve  pitch  settings. 

3.  We  coupled  the  3-D  uruteady  Navier-SuAes  solver  with  an  unsteady  full  potential 
flow  solver  to  arrive  at  a  hybrid  solver  which  solves  die  conqratationally  costly  Navier- 
Stokes  equations  only  in  die  small  viscous  region  surrounding  die  blade.  Elsewhere  die 
ineimnsive  potratial  flow  formulation  is  used.  The  formulation  is  fiiUy  duee-dhnensional, 
can  handle  unsteady  conqiressible  flows  with  embedded  shock  waves,  and  can  model  die 
effects  of  turbulence. 

4.  Many  of  the  2-D  and  3-D  analyses  have  been  modified  to  perform  aerodynamic 
design  of  ainoils  aiul  rotor  blades  as  well  We  tove  ma&  die  computer  codes  for  airfinL 
design  available  to  researchers  at  McDonh^jTSbnglaa'Helicopter  Co.  (Dy.  Ab"M»4 

'  Hassan),Befi  helicopter  Textron  (Jim  Narnuooie)  and  Sikatsl^OB^Moffitt). 

5 .  During  Jiis  period,  we  closely  worked  widi  odier  GttxpiL  Tech  researchers  on 
interdisciplin^  applications.  A  3-D  unsteady  panel  mediod,  mcotporaies  Dr.  Hodges' 
curved  beam  theory  was  developed  and  used  to  predict  the  leMl  lag  dan^g  characteristics 
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of  rotor  blades  in  hover.  We  worked  with  Dr.  Prasad  on  the  development  of  2-D  and  3-D 
unsteady  analyses  that  predict  the  effects  of  control  surface  motion  on  the  rotor  blade 
acoustics  and  vibratory  airloads. 


ERDcriments 

Team:  N.M.  Komerath,  S-G.  Liou,  and  Matt  Petrosld 

During  the  final  period,  Petrosld  wrote  up  results  for  submission  to  the  AIAA  Fluids 
conference,  and  prepared  a  paper  for  the  lichten  Award  Competition  placing  his  results  in 
tiie  context  of  our  previous  results,  CFD  results,  and  previous  work  on  both  helicopter 
blades  as  well  as  propellers  at  hi^  pitch  angles.  A  discussion  has  been  initiated  with 
researchers  at  NASA  Ames  (Peter  Talbot  and  (3iee  Tung)  on  die  general  observation  that 
propeller  blades  demonstrate  apparent  lift  coefficient  values  far  in  excess  of  their  predicted 
stalling  values  under  steady  rotating  conditicxis. 

Summary  of  fcyear-proiect 

The  experimental  ^rtion  of  this  task  set  out  to  study  the  flowfield  around  a 
complex  blade  tip  geometry  operating  at  high  pitch  angles.  This  was  an  experiment 
desired  to  assist  ctnnputational  research,  rxit  to  naodel  full-scale  phenranena.  Thus,  die  tip 
model  was  fabricated  using  the  surface  coordinates  used  in  the  Navier-Stokes 
computational  effort  A  single-blacted  rotor  was  used,  with  a  simple  inboard  geometry.  The 
operating  Reynolds  number  was  moderate,  within  reach  of  current  computational 
resources. 

Irutial  experiments  documented  die  behaidor  and  trajectory  of  the  tip  vortex.  The 
velocity  field  was  then  measured  using  IJDV,  at  low  and  hij^  pitch  angles.  The  results 
were  compared  directly  to  CFD  results  and  presented  in  a  1990  AHS  Forum  paper.  The 
rotor  experiments  were  dien  extended,  and  the  rotor  blade  was  then  used  as  a  fixed  wing  in 
the  wind  tunnel  at  a  freestream  speed  corresponding  to  the  tip  speed.  Velocity 
measurements  and  laser  sheet  visualization  were  then  u^  to  compare  flow  separation 
characteristics  and  spanwise  velocity  behavior  between  the  fixed-wing  and  rotor  cases. 
This  was  a  unique  experiment,  and  provided  results  that  supported  Ae  irxqjortance  of 
centrifugal  effects  in  the  rotor  flow  field.  The  results  were  presented  at  the  AIAA  Aerospace 
Sciences  Meeting  in  1991,  and  were  controversial.  Subsequently,  both  the  rotor  and  fixed- 
wing  experiments  were  extended  to  obtain  veloci^  field  with  finer  grid  resolution  near  die 
blatk  si^ace,  and  to  obtain  more  data  inboard.  Ihe  flow  visualization  was  also  extended 
using  a  new  smoke  generator.  This  providnl  a  surprising  benefit:  the  seeding  condensed 
and  deposited  on  the  rotating  blade  surface  in  the  regions  of  stagiiant  and  recirculating 
flow,  'ne  resulting  surface  flow  patterns  provided  excellent  correlation  with  velocity  field 
results  on  the  sinface  flow  behavior.  The  fixed  wing  experiments  were  also  repeated,  with 
the  fieestream  velocity  corresponding  to  the  tip  spe^  but  die  angle  of  attack  matched  to 
the  effective  rotor  angle  of  attack,  correcting  for  the  measured  inflow  velocity.  These 
experiments  provided  clear  visual  evidence  of  the  "inboard  notch  vortex"  that  appears  to 
control  flow  separatitm  on  the  tip  at  very  lugh  pitch  angles.  The  results  again  proved  that 
flow  separation  is  greatly  reducra  under  rotating  condititms  for  the  same  angle  of  attack. 
The  issue  of  centrifugal  pumping  remains  unresolved:  hypotheses  about  the  unambi^ous 
effects  of  centrifugal  punqiing  on  the  measurable  velocity  field  could  neither  be  ctmnrmed 
nor  rejected  based  on  the  measurements  obtained. 
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This  task  made  significant  contributions  to  our  measurement  and  visualization 
techniques.  The  copper  vapor  laser  was  made  operational  and  used  first  to  visualize  the 
vortex  structure  from  the  complex  blade  tip,  and  exanoine  the  reasons  for  the  observed 
unsteadiness  of  the  tip  vortex  trajectory.  The  initial  experiments  to  examine  the  concept  of 
Spatial  Correlation  Velocimetry  were  also  performed  here.  The  wax  deposition  method  of 
surface  flow  visualization  on  a  rotor  blade  was  observed  here,  and  is  expected  to  find 
extensive  use.  We  anticipate  strong  interest  in  the  problem  of  rotor  flow  at  high  pitch  in 
future.  The  "Himmelskamp  effect",  where  propellers  and  windmills  have  been  found  to 
operate  at  section  lift  coefficients  well  above  predicted  stall  under  steady  conditions,  is 
attracting  interest  at  NASA  Ames  (Peter  Talbot,  Cbee  Tung),  and  we  anticipate  increased 
collaboration  in  this  area. 
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II.  Aeroelasticity 


Task  i.  Unsteady  Aerodynamics  for  Rotor  Aeroelasticity 
D.  A.  Peters  and  D.  H.  Hodges 
Problem  Studied 

Rotor  aeroelastic  response  and  stability  analyses  have  traditionally  been  carried  out 
in  terms  of  modem  nonlinear  structural  dynamics  models  but  with  relatively  crude,  quasi¬ 
steady  aerodynamics.  For  the  hovering  flight  condition  this  combination  is  ^nerally  stated 
to  be  adequate,  although  without  proof.  In  forward  flight,  however,  it  is  generally 
recognized  that  this  type  of  formulation  is  severely  limited.  Incorporation  of  nonlinear, 
unsteady,  free-wake  analyses  based  on  panel  methods  into  such  codes  by  presently  known 
teclmiques  requires  the  introduction  of  hundreds  m  even  diousands  of  state  variables  and 
has  not  been  attenqited.  On  d^e  odier  hand,  relatively  sq^ust^aied  performance  codes'have 
existoi  for  some  time  that  make  use  of  duM-dimensi(m^  lifting  sunace  and  panel  methods 
for  determining  the  effects  of  tip  shape,  sweep,  etc.  on  thrust  and  performance.  However, 
such  codes  generally  do  ix)t  include  tte  effects  of  blade  aeroelastic  deformations.  It'is  clear 
from  existing  lifting  line  analyses  that  do  include  these  effects  that  they  could  have  a 
significant  inifluence  on  the  accuracy  of  calculated  danqring.  Alternatives  to  complete  ffee- 
wake  analysis  include  the  dynamic  inflow  theory  of  ^ters  and  co-workers  in  which  the 
induced  inflow  distribution  is  expanded  in  terms  (sf  an  arbitrary  number  of  radial  functions 
and  azimuthal  harmonics.  AnoAer  is  the  reduction  of  the  number  of  aerodynamic  state 
variables  in  the  analysis  by  expansicm  of  the  discretized  wake  parameters  in  terms  of  a 
relatively  small  number  of  functions  and  applications  of  some  variant  of  the  weighted 
residual  technique. 


Progress  During  the  Last  Reporting  Period 
a)  Hovering  StabUitv  and  Inflow 

Our  correlations  between  theoretical  danqring  predictitms  and  the  Army  dancing 
measurement  (Sharpe,  OTmiston,  et  al.)  have  be«n  gmtly  refined  and  improved.  First, 
we  have  revised  the  lift-curve  slope  and  drag  data  n^  by  previous  investigators  (and  by 
us)  through  the  use  of  NACA  tables  at  the  Reynolds  Nurnber  of  interest  (^0,000).  This 
change  has  improved  dancing  by  25%.  Second,  by  corrqraring  measured  redu^  flow  (as 
given  by  Sha^)  and  conqraring  with  predictions  from  nxrmentum  tiieory,  our  inflow 
model,  and  our  ptuiel  code,  we  have  determined  that  there  was  significant  recirculation  in 
the  test  chamber;  and  we  have  modified  our  results  to  reflect  tiiis.  Third,  detailed 
comparisons  Of  calculated  and  measured  fiequencies  have  indicated  that  the  torsional 
flexures  in  the  test  model  contributed  to  the  inplane  bending  flexibility.  Although  the 
percentage  change  is  not  large,  the  inplane  frequency  is  very  near  l.S/!rev;  and  inflow 
coiq)ling  is  significantly  affected  by  whether  ornot  the  finequency  is  greater  than  or  less 
than  1.5.  These  last  two  refinements  gave  anotiier  25%  in^aovement 

It  is  interesting  to  note  that  the  key  to  malting  us  kxA  more  deeply  into  static  inflow 
and  flexure  stiffness  came  frcnn  the  power  of  our  finite-state  nxxlel.  In  particular,  we 
noticed  that  the  tbeoreticd  collective-mode  danming  agreed  better  with  measured 
differential-mo^  daoqnng  than  did  the  theoretical  differential-mode  numbers.  We  tiien 


20 


Georgia  Institute  of  Technology 

Center  of  Excellence  for  Rotary  Wing  Aircraft  Technology 
Final  Report:  July  1, 1992  -  January  14,  1993 


looked  at  the  magnitude  and  phase  of  the  "lift  eigenvectors".  (Recall  that,  since  inflow  is 
represented  by  states,  eigenvector  information  contains  aerodynamic  as  well  as  structural 
information.)  The  information  in  the  lift  eigenvectors  led  us  to  believe  diat  small  changes  in 
frequency  about  1-SAev  as  well  as  small  changes  in  static  inflow  could  account  for  the 
damping  differences.  This  proved  to  be  true.  The  results  of  this  latest  correlation  will  be 
present^  at  the  49th  Annual  Forum  of  the  AHS.. 

b)  Forward  Flight 

Our  forward  flight  work  was  essential  completed  at  the  last  reporting  period.  We 
will  present  a  paper  on  the  results  at  the  coming  AIAA/AHS  DM  (inference  in  San  Diego, 
April  19-21. 

External  Interactions 

We  continue  to  work  with  Dr.  Robert  Qrmiston  of  the  Army  Aeroflightdynamics 
Directorate.  In  addition,  we  have  had  new  interactions  with  workers  at  Michigan  State  and 
with  the  Georgia  Tech  Flight  Simulation  Lab.  l^ally,  this  inflow  mc^el  is  being 
incorporated  in  ^QIAS-. 


Task  2.  Vibration  and  Trim  of  Elastic  Rotor  Blades  with  Dynamic 
Stall 

D.  A.  Peten 
Problem  Studied 

In  this  task,  we  have  been  studying  the  effect  of  dynamic  stall  on  rotor  trim,  on 
rotor  vibrations,  and  on  our  ability  to  compute  these  succes^ully.  Tliis  work  has  two  foci. 
The  first  is  on  development  of  good  anal^c  models  of  the  stall  mechanism.  In  this  area, 
we  began  with  the  ONERA  model  and  have  been  making  many  nxxlifications  to  it,  nx>st  of 
which  have  been  adopted  by  the  French  also.  The  second  focus  is  on  efficient  trim 
methodologies.  In  tlus  area,  we  have  worlmd  hard  on  harmonic  balance  techniques, 
periodic  shooting,  and  autopilots. 


Progress  Purine  the  Last  Reporting  Period 

We  have  re-visited  our  dynamic-stall  model  to  study  exactly  how  it  behaves  when 
coupled  with  our  finite-state  inflow  model.  Unfortunately,  however,  our  inflow  model  is 
3-D,  whereas  all  of  the  dyriamic  stall  data  are  2-D.  Thus,  we  first  had  to  develop  a  2-D 
version  of  our  finite-state  inflow  model.  This  was  done  under  a  grant  from  the  Army 
AerofLightdynarmcs  Directorate  (Robert  Ormiston,  Technical  Monitor).  We  have  now 
coupled  this  2-D  inflow  riuxlel  with  our  latest  version  of  dyruunic  stall,  and  this  is  giving 
excellent  data  conelation. 

In  trim  meAodologies,  we  have  yet  to  couple  successfully  periodic  shooting  and 
autopilot  in  a  hybrid  scheme.  However,  pure  shooting  and  pure  autopilot  are  both  mature 
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and  are  beiny  applied  in  aeroelastic  analyst  packages  at  NASA  Langley  Research  Center 
andliere  at  GERWAT/CXRT  foroptiiDizatum  qipUcadons. 

External  Interactions 

We  have  worked  with  NASA  Langley,  NASA  Ames,  the  University  of  Sao  Paulo, 
Advanced  Rotorcraft  Technology,  Inc.,  and  the  Washington  University  group.  Advance 
Rotorcraft  Technologies  is  putting  our  dynamic  stall  model  into  their  simuladtm  prop^mis. 


Task  3.  Unsteady  Aerodynamic  Testing  of  Model  Rotors 
G.  A.  Pierce,  N.  M.  Komerath,  and  S.G.  Liou 


Problem  Studied 

This  task  is  intended  to  identify  the  relationship  between  the  inflow  velocity  field  in 
the  vicinity  of  a  helicopter  rotor  blade  and  tire  associated  blade  pressure  distribution  for 
conditions  of  dynamic  pitch  control.  An  experimental  investigation  forms  the  basis  of 
determining  this  relationship  between  dynasac  inflow  and  unsteady  airloads.  Tests 
conducted  in  tire  AeroTech  facility  with  a  model  rotor  blade  consist  of  two  complementary 
investigations.  The  first  entails  the  measurement  of  absolute  pressure  distributions  on  the 
blade  surfaces,  while  the  second  involves  laser  velocimeter  measurement  of  the  inflow 
velocity  field  in  the  region  immediately  above  die  rotor  disk.  Recorded  synchronization  of 
these  measurements  with  respect  to  the  dynamic  pitch  ctmtrol  and  azimuthal  position 
permits  their  correlation  with  each  other  and  subsequent  identification  of  their 
interrelationship 


5-vear  Summary  of  Progress:  Experiments 

This  Task  sou^t  to  measure  tire  time-varying  inflow  to  a  rotor  excited  in  pitch,  and 
the  associated  effects  on  the  surface  pressure  field.  A  two-bladed  rotor  was  used  in  hover, 
with  pitch  excitations  driven  by  hydraulic  actuators  using  digitally-specified  4-per-rev 
harmonic  waveforms.  The  inflow  velocity  field  was  measured  using  LDV,  and  the  surface 
pressure  fluctuations  were  measured  using  surface-mounted  Kulite  absolute  pressure 
transducers.  Sevnal  majOT  problems  were  encountered  and  solved.  Some  of  these  were: 
acqmsition  of  velodty  ^ta  in  a  large  facility  with  high  enough  data  rate  to  capture  the 
detailed  waveforms  of  a  flowfield  excited  at  40Hz;  resolving  tiie  perturbations  due  to  a 
pitch  oscillation  of  1  degree  amplitude  in  a  periodic  flo^eld  with  large-amplitude 
fluctuations;  calibration  of  surface-mounted  absolute  transducers  to  enable  measurement  of 
extremely  small  pressure  changes;  transmission  of  data  from  extremely  low  pressure 
perturbations  from  tiie  rotating  fystem;  absolute  synchrcmization  of  the  systems  for  pitch 
excitation,  velocity  data  acquisition,  and,  later,  pessure  data  acquisition,  and  bter,  digital 
correction  of  the  pitch  excitation  waveform  for  the  complex  transfer  function  of  the 
excitation  systeno. 
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The  experimental  results  are  summarized  in  two  papers  submitted  to  the  AHS 
Journal;  a  third  pai^r  is  in  preparation  dealing  the  correlation  of  the  pressure  and  velocity 
results  with  analytical  methods.  In  addition,  the  velocity  field  results  have  been  used  by 
Ptof.  Peters  in  conelation  with  his  model  for  dynamic  inflow,  with  excellent  success.  This 
correlation  has  recently  appeared  iii  the  Journal  of  Aircraft.  Primary  conclusions  from  the 
velocity  field  experiments  are: 

Data  rates  over  4000  sec-1  have  been  achieved  in  backscatter  in  a  large  rotor  test 
facility.  The  scatter  in  individual  time  history  plots  is  minimal,  indicating  u^ormity  of 
seed  particle  size,  freedom  from  particle  lag  effects,  and  that  phase-averaging  is  not 
essential.  Agreement  of  the  steady-blade  data  with  analytical  results  demonstrates  that 
facility  flow  recirculation  errors  are  insi^ficant,  and  accurate  phase  synchronization  is 
achieved.  The  4  per  rev  inflow  variation  is  substantial,  even  for  1-deg.  oscillation 
amplitude.  The  variation  depends  on  excitation  phase,  and  is  relatively  insensitive  to  the 
blade  position.  The  inflow  response  to  pitch  excitation  decreases  in  amplitude  near  the 
blade  tip.  The  unsteady  inflow  variation  over  the  rotor  disk  was  successfully  synthesized 
assuming  axial  symmetry.  Substantial  phase  lag  effects  are  seen  in  the  flowfield.  The 
dynamic  inflow  variation  reveals  substantial  hysteresis,  implying  that  the  flow  is  totally 
tmsteady.  However,  the  unsteady  shedding  vorticity  can  treated  with  a  quasi-steady 
approach,  if  the  phase  lag  has  been  taken  into  consideration.  While  the  hysteresis  has  a 
radial  variation,  it  is  not  dnven  by  tip  effects. 

The  conclusions  finom  the  pressure  field  measurements  are: 

At  a  fixed  radial  and  chordwise  location,  a  180°  phase  shift  between  upper  and 
lower  surface  pressure  data  is  seen  even  at  low  amplitudes  of  pitch  oscillation. 
Instantaneous  perturbation  pressure  repeats  periodically  with  rotor  azimuth  angle,  as 
expected.  The  upper  and  lower  surface  penurbation  pressure  amplimde  profiles  remain 
close  to  each  other  and  converge  towards  the  trailing  edge.  Even  at  fixed  mean  pitch  angle, 
two  distinct  cases  exist  for  the  study  of  wake  dynamics,  depending  on  the  phase  between 
individual  blade  motions  and  thus  on  the  number  of  blades  and  n,  the  harmonics  of  the 
forcing  frequency.  At  a  fixed  radius,  amplitude  of  lifting  perturbation  pressure  increases 
with  the  reduced  frequenc>,  except  at  3-pCT-iev  and  S-per-rev  where  large  increase  occurs 
due  to  phase  effect  The  chordwise  distributions  of  the  phase  rotate  and  become  steeper 
with  reduced  fiequency.  The  reduced  frequency  decreases  as  we  go  outboard.  At  fixed  test 
conditions,  lifting  pressure  amplitude  deoeases  outboard.  The  phase  profile  rotates  and 
becomes  flatter.  The  above  two  effects  become  more  pronounced  with  increasing  reduced 
frequency.  At  high  reduced  frequency,  satisfaction  of  Kutta  condition  becomes  sluggish. 
Increasing  mean  pitch  increases  wake  smicing  and  affects  the  unsteady  pressure  amplitude 
near  the  leading  edge.  The  amplitude  ^t  increases  and  then  decreases  with  increasing 
mean  pitch  and  the  phase  angle  distribution  becomes  flatter.  At  the  blade  tip,  as  higher 
pitch  is  reached,  the  loss  is  sudden  and  severe. 
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III.  Structures 


Task  1:  Nonlinear  Beam  Theory 
D.  H.  Hodges 


Background; 

Many  existing  helicopter  blade  analyses,  including  that  of  the  Army’s  2GCHAS  program, 
are  based  on  simplifi^  beam  theories,  such  as  diat  of  Hodges  and  Dowell  (1974).  Inherent 
restrictions  in  such  analyses  that  are  inappropriate  for  realistic  rotor  blade  configurations 
include: 


(1) The  equations  are  restricted  to  moderate  rotation  theory  which,  for  certain 
problems  such  as  beaiingless  rotor  flexbeams,  is  inadequate. 

(2) They  do  not  treat  initial  curvature  and,  althou^  they  were  derived  to  treat  initial 
twist,  the  derivation  is  not  based  on  kinematics  appropriate  for  curvilinear 
coordinate  systems. 

(3) They  do  not  account  for  material  anisotropy,  shear  defonnation,  or  warping,  all 
of  which  are  essential  for  oxxipositB  blade  moling. 

(4) They  are  based  on  a  Green  strain  formulation  which  produces  man^f  sup^uous 
terms  in  the  equations  of  motion.  Based  on  past  ad  hoc  methods  of  sinqilification 
(such  as  ordering  schemes),  removal  of  these  terms  is  not  straightforward  nor  does 
it  result  in  practical  sixx^tliBcations,  since  die  resultin,,  equations  fill  several  pages. 

Originally  we  set  out  to  overcome  these  problems,  and  we  believe  that  we  have  achieved 
this  and  more.  Some  aspects  of  our  theory  are  quite  mature  now,  and  others  are  in  die  final 
stages  of  being  developed.  Still  odiers  are  hcyaid  the  scope  of  tte  present  project  and  must 
be  developed  under  the  auspices  of  separately  funded  projects.  Also,  further  validation 
needs  to  be  done,  much  of  which  remains  in  planning  stages  because  earlier  graduate 
research  assistants  were  lost  through  the  PhJ).  qualifying  exam. 

Summary  of  Theorv  Developed  to  Date: 

The  work  accomplished  to  date  has  led  to  a  modular  theory  which  we  view  as  ideal  for 
rotor  blade  modeling  as  well  as  for  other  plications.  Our  theosy  is  based  cm  a  decoupling 
of  the  sectional  andysis  (to  obtain  elastic  cross-sectional  properties)  from  the  globd 
deformation  analysis. 

Berdichevsky  (1981)  appears  to  be  the  first  in  the  literature  to  show  that  consistent 
application  cd  die  variational-asyn^totical  method  to  die  ^uations  of  three-ddmensional 
geometrically  nonlinear  elasticity  for  beams  leads  to  two  silkier  problems:  a  linear,  two- 
dimensional  cross-sectional  analysis,  from  which  the  sectional  prcqierties  can  be  calculated, 
and  a  nonlinear,  one-eHmensional  beam  analysis  which  uses  those  properties  in  calculation 
of  the  beam  deflections,  rotations,  etc.  The  sectional  analysis  also  yields  a  set  ttf  influence 
functions  which  can  be  used  to  calculate  the  three-^mensional  stress,  strain,  and 
displacement  throughout  the  beam  once  the  cme-dimenskmal  problem  is  sdved 
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Accordingly,  we  have  constructed  our  theory  in  a  modular  fashion  with  four  parts: 

(a)  exocf,  nonlinear,  Idnematical  equations  for  beam  displacements  and  rotations  in 
terms  of  intrinsic  strain  measures  (stress  resultants  do  not  appear  in  these 
equations); 

(b)  exact,  nonlinear,  intrinsic,  global  equilibrium  equations  in  terms  of  physical 
stress  resultants  and  intrinsic  strain  measures  (displacement  and  rotation  variables 
do  not  appeal  in  these  equations); 

(c)  asymptotically  correct  constitutive  equations  which  relate  the  physical  beam 
stress  resultants  (section  force  and  moment  measure  numbers)  to  the  intrinsic  strain 
measures  (displacement  and  rotation  variables  do  not  tqipear  in  these  equations); 

(d)  influence  fractions  which  give  asymptotically  correct  expressions  for  all  three- 
dimensional  displacement,  strain,  and  stress  variables  in  terms  of  one-dimensional 
intrinsic  strain  measures  (one-dimensional  stress  resultants,  displacement  and 
rotation  variables  do  not  appear  in  these  relations). 

The  first  three  modules  of  the  thec^  can  be  combined  to  form  a  single  mixed  variational 
statement  in  space-time,  described  in  detail  by  Hod^s  iintemational  Journal  cf  Solids  and 
Structures,  1990).  All  four  parts  must  account  for  initial  curvature  and  twist  In-  and  out- 
of-plane  St-Venant  warping  is  a  by-product  of  (c)  and  is  used  to  find  (d);  however,  no 
Wiping  variables  appear  in  any  of  the  final  equations.  Parts  (a)  and  (b)  along  with  the 
resitits  of  (c)  constimte  a  closed  theory  which  can  be  put  into  a  very  compact  form  witirout 
the  necessity  of  ordering  schemes  or  other  ai^ximations.  Part  (d)  is  used  for  obtaining 
the  complete  displacement  field  and  die  strain  or  stress  field  throughout  the  beam  once  the 
equations  in  (a)  —  (c)  are  solved.  Certain  aspects  of  parts  (c)  and  (d)  are  complete,  but 
some  are  still  under  develq)ment 

It  should  be  noted  that  treatment  of  the  restrained  warping  requires  additional  variables  of 
all  types.  This  is  very  much  analogous  to  nxxleling  a  flexible  body  with  rigid-body  modes 
and  flexible  modes,  the  St-Venant  solutions  being  analogous  to  the  rigid  body  mo^s.  The 
edge-zone  behavior  is  governed  by  higher-order  deformation  modes  of  the  section,  very 
much  analogous  to  higher-frequency  modes  in  deformable  bodies.  The  treatment  df  edge- 
zone  phenomena  and  high-finquency  d^amics  must  be  regarded  as  coupled,  as  previous 
work  has  shown.  Primarily  because  of  its  con^lexity,  further  development  on  this  aspect 
of  the  problem  is  outside  tiie  scope  of  the  present  jroject  A  separate  prc^osal  for  research 
on  this  subject  has  been  written  by  Dr.  V.  L.  Berdichevsky  (one  of  the  leading  authorities 
in  the  world  on  this  subject)  and  the  P  J.  and  submitted  to  die  ARO. 

Constitutive  Equations  and  Influence  Functions  fCesnikV 

Motivation:  All  the  present  work  is  being  devoted  to  develqnng  a  practical,  con^tationally 
efficient  method  to  extract  the  elastic  cross-sectional  properties  of  an  anisotropic  beam  wim 
arbitrary  cross-sectional  materials  and  geometry.  The  most  systematic  framework  to  our 
knowledge  is  the  variational-asymptotical  method  of  Bertfichevsky.  This  analysis 
methodology  is  able  to  provide  the  asynqitotically  contet  influence  functiras,  allowing  us 
to  express  the  three-dimensional  field  variables  in  terms  of  ons-dimensional  generalized 
strain  measures.  This  means  that  for  arbitrary  nonlinear  deformation  of  anisotropic  beams, 
we  ^  able  to  determine  the  displacement,  strain,  and  stress  throughout  the  beam. 
Published  work  prior  to  ours  which  applies  tins  method  to  beams  does  not  consider  the 
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generally  anisotropic  beam  beyond  the  first  approximation  (clarsical  theory  for  bending, 
torsion,  and  extension),  nor  are  these  method  developed  into  computational  methods 
developed  to  calculate  the  elastic  constants.  Indeed,  Aeie  is  a  long  way  between  the 
mathematical  theory  and  a  practical  computational  methodology  for  realisdc  rotor  blade 
structures. 

Accomplishments:  We  developed  the  equations  for  the  first  approximation  for  anisotropic 
beams  section  constants  and  programmed  for  solution  by  using  the  finite  element  method  to 
disciedze  the  beam  cross  secdoiL  A  finite  element  code  was  written  to  evaluate  the  warping 
and  the  stiffness  matrix  for  a  beam  cross  secdon  comprised  of  laminated  orthotropic 
materials  in  which  the  fiber  direcdon  is  arbitrary.  Note  thM  our  analysis  is  not  restricted  to 
beams  which  have  sections  of  any  particular  geometry  (such  as  thin-walled).  The  element 
used  is  a  four-noded  planar  rectangular  element  with  three  degrees  of  fioedom  per  node. 
Element  matrices  are  formulated  by  exact  integradon  using  Mathematica  (a  symbolic 
manipulation  program).  The  code  is  quite  portable  and  can  t«  run  on  deslctop  computers 
such  as  PC’s  and  Macintoshes.  It  is  refened  to  as  VABS  (Variadonal-Asymptodcal  Beam 
Sectional)  Analysis. 

Progress  During  This  Reporting  Period:  In  this  reporting  period  we  have  extended  the 
previous  work  to  include  die  first-order  effects  of  initial  curvature  and  twist.  The  warping 
m  this  case  has  two  parts:  the  first  is  the  same  as  that  for  die  prismatic  case  and  the  second 
is  a  first-order  correction  to  account  for  initial  curvature  and  twist  The  corrected  warping 
can  be  used  to  find  accurate  three-dimensional  displacement  and  strain  fields  within  the 
beam,  but  the  correction  of  the  warping  is  not  needed  to  find  the  elastic  constants.  The 
theory  is  such  that  the  only  static  coupling  between  these  new  degrees  of  freedom  and  the 
existing  six  is  elastic,  i.e.,  in  the  constitutive  law.  The  theory  for  initial  twist  and  curvature 
was  reported  in  papers  at  SECTAM  and  DM.  Numerical  results  obtained  to  date  show  a 
quadratic  convergence  rate  against  the  number  of  elements. 

Future  Work:  Our  present  element  is  not  sufficient  for  realistic  blade  noodeling,  and  plans 
are  being  made  to  expand  the  element  library  to  include  isoparametric  elements,  which  treat 
curved  r^oil  surfaces  more  accurately  than  the  present  method  can.  We  also  are  ma^g  an 
effort  to  improve  the  computational  efficiency  of  the  method.  There  are  several  instances  in 
which  certain  mathematical  transformations  which  are  appropriate  for  analysis  (to  obtain  a 
well-posed  minimum  principle,  for  example)  are  not  necessarfiy  the  best  ones  to  use  in 
computational  algorithm. 

In  the  process  of  extending  our  theory  and  code  to  account  for  initial  curvature  and  twist,  a 
great  deal  of  progress  was  made  in  ^veloping  a  frameworic  for  incorporating  degrees  of 
fieedom  other  tlum  the  six  rigid-body  translational  and  rotational  degrees  of  fieedom  for  the 
cross  section  now  treated.  This  would  be  important  for  situations  in  wltich  the  warping 
must  be  restrained,  such  as  in  bearingless  rotors.  The  theory  is  such  that  the  only  static 
coupling  between  Aese  new  degrees  of  freedom  and  the  existing  six  is  elastic,  i.e.,  in  the 
constimtive  law. 

Vahdation  (ReeceV 

Another  aspect  of  our  research  calls  for  the  validation  of  the  cross-sectional  analysis.  We 
have  now  compared  our  numerical  results  with  data  fix>m  composite  beam  experiments 
done  at  Massachusetts  Institute  of  Technology  and  at  the  University  of  Maryland  C)ur  own 
experimental  work  is  progressing  slowly.  Rather  than  duplicate  what  oAers  have  done,  tiie 
focus  of  our  work  will  be  to  address  whether  spanwise  nonumformity  significantly  affects 
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the  elastic  constants,  starting  with  the  isotropic  case.  Our  first  experimental  results  should 
be  available  during  the  next  reporting  period,  if  all  goes  according  to  present  plans. 

Numerical  Work  rChune  and  ShaneV 

In  order  to  apply  (a),  (b),  and  (c)  efficiently,  it  is  necessary  to  take  advantage  of  the 
sparsity  of  the  coefficient  matrices  and  to  try  to  optimize  sparsity  by  choice  of  shape 
functions.  To  date  all  our  nonlinear  finite  element  work  has  centered  on  the  use  of  the 
weakest  form  of  the  equations  and  the  crudest  possible  shape  functions.  In  this  approach 
numerical  element  quadrature  is  not  necesst^  if  the  properties  and  forcing  functions  are 
integrable  functions.  However,  for  non-uniform  beams  we  do  not  know  whether  it  is 
necessary  to  integrate  the  properties  exactly  or  if  it  would  be  ^st  as  accurate  to  approximate 
them  as  piece  wise  constant  to  go  along  with  the  sha^  functions.  This  has  a  dirrct  bearing 
on  how  one  approaches  development  of  p-version  finite  elements  in  the  nonlinear  case.  We 
know  that  these  elements  perform  very  well  relative  to  the  crude  approach  if  the  element 
properties  are  constant,  the  forcing  function  is  constant,  and  the  problem  being  treated  is 
linear.  Although  what  happens  in  the  general  case  is  still  being  studied,  we  have  now 
established  that  higher-orda  elements  do  inqjtove  the  accuracy  si^^cantly  for  beams  with 
constant  and  linearly  varying  properties.  The  next  step  is  the  examine  what  happens  when 
numerical  quadrature  is  us^  instead  of  exact,  and  to  go  through  the  same  jvocess  for  a 
nonlinear  problem.  Further  unansweml  questitms  residb  in  the  treatment  of  time-domain 
issues.  How  does  one  ensure  numerical  stability  in  time-marching  algorithms  for  mixed 
models  of  beam  dynamic  behavior?  We  view  this  as  outside  the  scope  ctf  the  current  projea 
and  plan  to  address  it  later  under  the  auspices  of  a  new  project 

Related  Research 

Oir  work  has  been  applied  to  the  dynamics  and  aeroelastic  stabili^  analysis  composite 
rotor  blades  in  hover  under  odier  ARO-sponsorship;  this  activity  is  in  a  no-cost  extension 
and  will  end  in  September  1992.  Our  model  appears  to  be  well-suited  for  low-fiequency 
aeroelastic  stability  analysis.  Our  work  to  extend  apply  similar  methodology  to  laminated 
plate  behavior  is  being  funded  by  the  Army  Aaostructures  Directorate  at  NASA  Langley 
Research  Center.  Our  being  able  to  work  these  problems  in  parallel  has  been  extremely 
beneficial. 

External  Interactions: 

with  Drs.  Robert  A.  Ormiston  and  Gene  C.  Ruzicka  of  Aeroflightdynamics  Directorate 
concerning  GRASP  and  2GCHAS  (PI  spent  three  weeks  at  AFDD  during  Summer  1991) 

with  Dr.  Charles  Rogers,  Dr.  Gene  Sadler  and  Mr.  Mark  Dreier  of  Bell  Helicopter 
concerning  finite  element  modeling  of  rotor  blades 

with  Prof.  Larry  Rehfield,  University  of  C^alifomia,  Davis  concerning  extension  of  his 
sectional  analysis  theories  and  their  qiplication  to  ^ledfic  configurations 

with  Dr.  T.  Kevin  O’Brien  and  Mr.  Howard  Hinnant  of  the  Army  Aerostructures 
Directorate  and  Dr.  Alexander  Tessler  of  NASA  Langley  concerning  laminated  plate 
problems 


Task  2.  Deleted 
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Task  3.  Rotorcraft  Vibrations  and  Structural  Dynamics 
S.  Hanagud  and  J.  L  Craig 
Problems  Studied 

As  a  part  of  this  task,  we  conducted  research  in  two  different  but  related  areas  of 
structural  dynamics.  The  subject  matter  of  the  first  part  was  smart  or  adaptive  structures. 
The  subject  matter  of  the  second  part  was  structural  dynamic  system  identificatitni. 

SMART  STRUCTURES 

Smart  structures  can  be  defined  as  a  class  of  structures  that  have  built-in  adaptive 
capability  and/or  intelligence.  Such  an  adaptive  ct^abUity  or  intelligence  is  used  to  optimize 
die  performance  of  a  structural  system  to  changing  environments.  The  subject  matter  is  an 
intodisciplinary  field.  In  our  experience,  the  research  has  involved  a  combination  of  the 
fields  of  structural  mechanics,  control,  fluid  mechanics,  and  mathematics.  The  research 
has  been  both  theoretical  and  experimental. 

Vibration  Control 

Our  work  started  during  the  second  year  of  the  first  phase  of  CERWAT.  We 
initiated  the  work  on  mathematical  modeling  and  optimum  vibration  control  of  beams  with 
bonded  piezoceramic  sensors  and  actuators.  We  followed  this  work  widi  smart  structural 
designs  that  can  provide  robust  control.  The  robustness  was  provided  to  account  for  some 
debonding  of  sensors  and  delaminations  that  would  alter  the  stiffness  and  differential 
equations  of  the  model.  We  developed  the  needed  equatitms.  Controllers  were  designed 
by  using  H-infinity  and  Mu-synthesis  techniques.  We  have  also  ^veltmed  techniques  to 
control  nonlinear  vibrations.  A  first  phase  of  e;q)eriments  to  control  noitimear  vibratitms  is 
being  develqped. 

Smart  Structures  in  Higher  Hannonic  Control 

We  have  completed  a  feasibility  study  in  this  area.  The  objective  of  the  feasibiliQr 
study  was  to  reduce  vibrations  at  select^  locations  of  an  aiifirame  of  a  rottneraft  by  the  use 
of  smart  structures  concepts.  The  usual  higher  harmonic  control  considers  a  reduction  of 
n-per-rev  oscillator  loads  transmitted  to  the  airframe  fix>m  the  rotor  systenL  In  our  woik, 
we  have  used  distributed  sensors  and  actuators  in  the  airfranoe  and  an  H-infinity  controller 
to  control  higher  harmonic  components  of  vilnations  at  selected  locations  of  the  airframe. 
We  are  also  studying  a  combination  of  rotor  mounted  systems  and  airframe  mounted 
sensor/actuaior  system  to  optimally  reduce  vibrations  at  selected  locations  of  a  helicopter. 

Smart  Structures  in  the  Health  Monitoring  of  Structures 

An  mqroi^t  application  of  smart  structures  is  in  tiie  area  of  healtii  noonitcaring  of 
structures.  A  typical  health  monitoring  system  should  be  able  to  identify  Ae  damage  when 
it  occurs  and  tten  automatically  control  the  growth  of  dunage.  We  have  completed  a 
feasibility  study  of  delamination  detection  arc  control  of  delamination  growth,  fo  thl . 
feasibilify  study,  Ae  delamination  detectitm  techniqw  is  based  on  tiie  observed  change  in 
Ae  dynamic  response  of  a  delaminated  structure.  The  automatic  control  of  die  growA  of 
Ae  delanodnaticHis  is  accmiqilished  by  reducing  the  magnitude  of  interlaminar  stresses  by 
using  Aeconcqpt  of  smart  structures  and  a  control  moment  feedback.  Our  next  step  will  be 
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the  validation  of  our  work.  We  have  experimentally  completed  a  first  phase  of  the 
validation  of  delamination  detection. 

Adaptive  Airfoils  and  Smart  Structures 

The  ability  to  actively  change  the  shape  of  an  airfoil  section,  while  it  is  in  motion, 
has  been  the  dream  of  many  aerodynamicists,  dynamicists,  and  control  engineers  for  years. 

At  Georgia  Tech,  we  have  developed  four  concqits  to  inake  such  active  changes  in  the 
airfoil  shape.  Our  first  attempts  were  based  on  the  use  of  shape  menxny  alloys.  We  were 
able  to  get  large  leading  and  trailing  edge  deflections  by  the  use  of  shape  memory  alloys. 
However,  because  of  the  needed  heating  and  cooling  cycles,  we  had  difficulties  in 
obtaining  the  needed  frequency  responses.  This  difficulty  has  been  resolved  and  we  are 
now  in  the  process  of  developing  a  two  bladed  rotor  system. 

As  a  next  step,  we  have  used  piezoelectric  stacks  to  obtain  flap-motion  in  the  range 
of  3-10  Hz.  This  is  based  on  a  canttiever  beam  actuation.  Angles  in  the  range  of  2  to  S 
degrees  have  been  achieved.  We  are  now  studying  methods  of  using  piezoelectric  actuators 
to  achieve  higher  frequency  ranges  and  larger  angular  deflections  without  the  use  of  the  , 
cantilever  beam  concept 

We  have  also  theoretically  studied  tite  use  of  adaptive  airfoils  in  actively  controlling 
the  blade-vortex  interaction.  During  tiie  operation  ci  a  rotorcraft,  rotOT  blades  interact  wiA 
vortices  shed  by  preceding  blades.  As  a  result  of  the  interaction,  lar^  pressure  pulses  are 
created  at  the  leading  edges  of  the  airfoil.  We  have  shown  the  feasibility  of  rejlucing  the 
magnitude  of  these  large  pressure  pulses  by  tiie  use  of  adaptive  airfoils.  We  have  used 
smart  actuation  to  change  the  camber  of  the  airfoil  and  optimum  control  techniques  to 
reduce  die  pressure  magnitudes.  Work  is  ctmtinuing  in  this  area  under  a  SBIR  with  the 
U.S.  Army  AFDD  at  Atks  Research  Center. 

Jitter  Vibrations  in  Rotating  and  Pointing  Systems 

The  concept  of  smart  stroctures  has  also  been  implemented  for  the  jitter  vibration 
reduction  of  flexible  pointing  and  rotating  systems.  Tracking  a  desired  trajectory  using 
flexible  rotating  and  pointing  systems  is  a  nonlinear  control  problem  due  to  the  combined 
effect  of  large  rigid  body  motions  and  small  flexible  body  motions.  Smart  structures 
concepts  have  been  used  in  the  form  0(  distributed  sensors  and  actuators  to  reduce  jitter 
vibrati^ons  caused  by  the  flexible  body  motions.  Hanagud,  et  al.  used  piezoceramic 
sensors  and  actuators  for  this  jitter  reduction  ^licatioiL 

Micmmotora  as  Smart  Actuators 

Actuators  of  miniature  size  are  called  micromotors.  Micromotors  are  either 
electrostatic  motors  or  motors  developed  using  smart  materials  like  piezoelectric  materials, 
electrostrictive  materials,  shape  niemoiy  alloys,  etc.  One  of  the  motivations  for  the 
develcpment  dT  nucromotors  is  based  on  the  fact  that  smaller  systems  can  move  small  parts 
much  mter  than  larger  systems.  Moreover,  thermal  expansion  and  vibration  problems  are 
minimal  when  smaller  systems  are  used.  In  addition,  it  is  easier  to  obtain  accuracy 
with  smaller  systems  and  because  the  flom*  ipace  required  is  minimal,  mey  can  be 
incorporated  at  moicate  locations  widiout  much  difficult. 

At  Georgia  Tech,  AUen  and  his  colleagues  have  developed  a  novel  method  of 
fabricating  high  aspect  ratio  electrostatic  micromotors.  This  process  used  polyimide 
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electroplating  forms  for  the  fabrication  along  with  standard  commercially  available 
materids.  The  fabrication  process  begins  with  the  making  of  an  oxidized  silicon  wafer 
substrate  coated  with  a  four  layer  metal  system.  Chromium  is  used  as  the  first  metal  and 
acts  as  an  adhesion  layer.  Copper  is  used  for  the  second  metal  for  plating  die  stator  and  pin 
components  of  the  motor.  Again,  a  layer  of  chromium  is  used  for  metal  three  for  the 
purpose  of  protecting  the  copper  coating.  This  layer  of  chrcnnium  will  also  act  as  a  release 
layer  for  the  rotor  smicture.  A  layer  of  copper  us^  for  fhe  fourth  metal  acts  as  a  seal  layer 
for  the  rotor  structures.  We  have  completed  a  preliminary  study  to  incorporate 
micromotors  as  smart  actuators. 

IDENTIFICATION 

The  field  of  structural  djmamic  system  identification  has  been  an  active  field  for  the 
past  two  decades.  However,  with  the  exception  of  some  preliminary  work  by  Meirovitch 
and  his  colleagues,  there  has  been  v^  litde  effort  devoted  to  development  of  iden^cation 
techniques  that  can  relate  an  identified  mathematical  mo^l  to  changes  in  the  physical 
variables  and  can  apply  the  method  to  large  scale  helicopter  airframe  structural  models 
(AH-IG  airframe  finite  element  models).  A  part  of  the  support  for  this  work  was  also 
provided  by  NASA.  Our  method  is  a  two*step  method  and  can  be  designed  to  .consider 
specific  variables  like  damping  constants,  dynamic  stiffiiesses  or  boundary  conditions. 

A  second  area  that  is  of  significant  interest  in  helicopter  rotor  dynamics  and 
structural  dynamics  is  the  area  of  nonlinear  system  identification.  As  a  part  of  the  first 
phase  of  the  CERWAT  program,  we  devel^d  a  perturbation  method  for  nonlinear 
structural  dynamic  system  identification.  During  the  second  phase,  we  have  develop^ 
non-linear  structural  dynamic  system  identification  techniques  on  the  basis  of  Hammerstein 
integral  operators.  Our  method  is  call  the  Hamnterstein's  feedback  model  for  nonlinear 
system  identification.  This  method  has  resulted  in  practical  techniques  tiiat  can  be  applied 
to  multi-degree-of-freedom  systems. 

TECHNOLOGY  TRANSFER 

Our  program  has  included  technology  transfer  with  the  Army  laboratories  at  Fort 
Eustis,  NASA-Langley  Research  Center,  and  NASA-Ames  Research  Center.  We  have 
also  worked  with  SikOTsky  aircraft  in  the  areas  of  structural  integrity  and  the  development 
of  smart  airfoils.  We  have  worked  with  Bell  Helict^ter  in  the  area  of  structural  dynamic 
system  identification. 

As  a  pan  of  the  technology  transfer  program,  we  have  published  27  papers  in  the 
area  of  our  research. 
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Task  4.  Damage  Resistance  in  Rotorcraft  Structures 
E.A.  Annanios 


Problem  Studied 

The  piimaiy  objective  of  this  research  is  the  develcqmient  of  a  damage  resistance 
ctmcept  in  cmnposite  rotorcraft  and  airframe  structures  by  tailoring  the  microstructure. 

Final  Report: 

This  repon  is  an  overview  of  the  research  performed  in  this  task.  It  describes  the 
major  accomplishments,  the  lessons  learned  and  die  future  investigation  which  stems  from 
its  findings. 

Damage  resistance  through  microstructural  tafloring: 

This  concept  is  based  on  the  fact  that  a  given  damage  mode  in  laminated  composite 
does  not  occur  in  isolation  but  is  often  accon^anied  by  other  noodes.  Microstructural 
tailoring  takes  advantage  of  the  interaction  a[  the  different  damage  modes  to  create  a 
resistance  behavior.  The  theoretical  basis  of  this  conc^t  was  developed  in  Ref.  1  and  its 
validation  was  provided  by  testing  ply  drop  ctmfiguratitxis  which  were  designed  based  on 
the  microstructure  tailoring  concept  The  design  was  intended  to  create  a  resistance  to 
delamination  or  ply  separation,  a  pnma^  damage  mode  in  laminated  coo^posite.  Test  have 
verified  delamination  resistance  behavior.  The  average  load  at  failure  relative  to  the  onset 
of  delamination  is  1.5.  This  corresponds  to  a  relative  fracture  touglmess  increase  of  1.22. 
Delamination  resistance  behavior  was  created  by  interaction  of  matrix  microciacking  with 
delamination.  Evidence  of  matrix  microcracking  with  delamination  growth  was  provided 
by  scanning  electron  miciography  and  synchrotron  microradio|T^>hy  [2]. 

The  inclusion  of  an  a^esive  la^r  or  a  resin  rich  pocket  results  also  in  a 
delamination  resistance  behavior.  A  quantitative  fiactogrqthy  noodel  [5]  was  developed  in 
order  to  explain  the  observed  delammatitm  arrest  pheixnmia.  The  variations  in  fracture 
surface  morphology,  found  by  the  quantitative  fia^graphy  and  statistical  analysis,  were 
determined  to  be  caused  by  variations  in  the  thickness  of  the  interlaminar  matrix  layer, 
rather  than  by  events  during  fracture.  The  laminates  were  unidirectional  ply  drop  mada  of 
brittle  as  well  as  toughenra  material  systems.  •  A  simple  model  has  bwn  proposed  that 
explains  the  arrest  phenomemm  in  terms  of  one-dimensional  model 

Fracture  Analysis  of  Transverse  Crack-tip  and  Ftee-edye  Delamination; 

Since  the  damage  resistance  concept  uses  the  interaction  of  damage  modes,  an 
investigation  of  transverse  crack-t^  and  free  edge  delamination  was  performed  [4].  Critical 
'  loads  and  delaminatitm  modes  were  identified  and  compned  with  experimental  results. 
Hygrothermal  effects  were  included  to  make  die  comparisons  realistic.  Hygrodiermal 
stresses  due  to  the  cure  cycle  can  have  a  significant  infiuenoe  on  the  delaminuxm  behavior 
of  laminated  composites.  This  was  shown  in  Ref.  5  for  die  interiaminar  stresses  and  strain 
energy  release  rate  associated  widi  Mode  I  fiee-ed«  delamination. 

The  investigation  of  damage  modes  in  laminated  conqxjsites  has  focused  on 
extensitm  loading,  however,  rotororaft  structures  are  subjected  to  combined  extension, 
bending  and  torsion  loads.  In  order  to  ensure  damage  tolerance  of  rotoraaft  structures 
combined  loading  effects  should  be  ccmsidered.  An  inteilanmnar  fracture  analysis  for 
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laminates  subjected  to  combined  extension,  bending  and  torsion  loads  has  been  developed. 
This  analysis  is  found  to  be  also  applicable  to  unsymmetiical  laminates. 

Analysis  of  Laminated  Composites  under  P/^mbined  Loading 

An  analytical  model  for  laminates  subjected  to  torsion  was  developed  for 
unidirectional  and  cross-ply  laminates[6].  Tbis  anal^s  was  extended  to  include  combined 
bending  and  torsion  in  Ref.  7  and  combined  extension,  bending  and  torsion  loads  in  Ref. 
8. 

It  is  found  that  this  analytical  metiiod  can  also  be  applied  to  unsymmetiical  laminate 
under  uniform  extension.  This  is  of  significant  practical  implications  as  damage  modes 
often  alter  the  initial  symmetry  existing  in  laminate  designs.  Moreover,  unsymmetrical 
laminate  designs  are  also  used  for  elastic  tailoring.  Potential  delamination  sites  in  a  [-  q  / 
(90  -  q)2  /  -q/q/(q  -  90)2  /  q  iT  class  of  laminates  were  predicted  based  on  the 
interlaminar  stresses  and  the  tohd  strain  energy  release  rate.  This  class  of  laminates  is 
designed  to  exhibit  extension-twist  coupling  wi&  no  initial  waiping  due  to  curing  stresses. 

A  computational  scheme  was  developed  in  order  to  predict  the  onset  of 
delamination.  The  analysis  involves  two  steps.  In  the  first  step,  potential  delamination 
sites  or  critical  interfaces  are  determined.  This  is  done  by  assuming  delaminations  at  the 
free  edges  of  one  interface  at  a  time  and  predicting  the  sign  of  the  interlaminar  p^l  stress  in 
the  neighborhood  of  tiie  delamination-tip.  Tl»  closed  form  expressions  i^vided  by  the 
analysis  makes  it  possible  to  perform  this  procedure  at  minimum  computational  effort  A 
positive  peel  stress  represents  a  potential  critical  interface.  In  the  second  step  a  complete 
analysis  of  the  laminate  is  performed  with  delamination  at  the  free  edges  of  the  cntical 
interfaces.  The  predicted  delamination  sites  based  on  the  interlaminar  peel  stress  is 
consistent  with  the  total  strain  energy  release  rate  prediction. 

The  results  of  this  research  are  presented  in  Ref.  9 .  To  the  best  of  our  knowledge 
this  is  the  first  closed  form  fracture  analysis  model  fcr  unsymmetiical  laminates. 

The  potential  benefits  of  applying  extension-twist  coupled  laminates  in  rotorcrafr 
composite  structures  provided  tiie  motivation  for  assessing  the  [-  q  /  (90  -  q)2  /  -q/q/(q’ 
90)2  /  q  ]T  class  of  laminates  and  developing  means  for  maximizing  coupling  under 
mechanical  loads  while  ensuring  hygrotheimal  staUlity. 

Analysis  of  Laminates  with  Optimum  Extension-Twist  Coupling 

The  shear  defonnation  model  devdoped  for  unsympietiical  laminates  was  qiplied  to 
the  desi^  and  aiudysis  of  laminates  with  optimum  extensitm-twist  coupling.  This  is  also  a 
prer^uisite  to  assessing  the  effect  of  damage  in  laminated  composites  designed  for 
maximum  coupling.  A  constrained  optimization  scheme  was  developed  in  order  to 
determim  die  staci^g  s^tienoe  that  marimiges  the  extension-twist  coupling  in  a  composite 
laminate  while  maintaining  its  hygrotheimal  stability  [10\.  The  results  were  conqnm  witii 
the  hygrothermally  stable  laminates  obtained  by  stacking  a  set  of  rotated  [0/90]  plies  such 
asthe[-q/(90-q)2/'q/q/(q  -  90)2  /q  ]T  class  of  laminates  .  It  was  found  that  the 
laminate  coi^gurations  resulting  from  the  optimization  scheme  have  inqtroved  coupling. 
Their  coupling  magnitude  however,  is  voy  sensitive  to  changes  in  fiber  ang^e.  This 
restricts  their  implementation  in  practical  designs. 

Eight  sets  of  laminates  with  varying  extension-twist  coupling  were  made  in  house 
with  Hercules  AS4/3502  Graphlte/Epoxy  and  tested  using  a  custom  loading  transducer. 
The  laminates  test  results  were  compared  with  analytical  predictions  horn  tite  analytical 
model  developed  and  firom  Classicil  Lamination  Theory  ((XT).  It  was  found  that  shear 
deformation  has  a  neglipble  influence  oa  the  extension  twist  coupling  and  predictions  from 
both  theory  were  in  close  agreement. 


32 


Georgia  Institute  of  Teclmology 

Center  of  Excellence  for  Rotary  Wing  Aircraft  Technology 
Final  Report  July  1, 1992  -  January  14,  1993 


The  result  of  this  work  indicate  diat  die  hygrothennally  stable  laminates  obtained  by 
stacking  a  rotated  0/90  group  of  plies  provide  optimum  extension-twist  coupling  and 
robustness.  Moreover,  the  extension-twist  coupling  in  this  class  of  lanoinates  is  accurately 
predicted  by  Qassical  Lamination  Theory. 

Effect  of  Free-edee  Delamination  on  Extension-twist  Coupling: 

The  optimum  configurations  achieved  in  Ref.  10  are  the  basis  for  assessing  the 
influence  of  damage  on  elastically  tailored  laminates.  Eight  sets  of  laminates  with 
embedded  free-edge  damage  were  fabricated  in  house  with  Hercules  AS4/3S02 
Graphite/Epoxy  and  [30/-6(]^0/-30/602/-30]t  layup.  The  damage  was  simulated  by 
placing  a  thin  Teflon  I%P  film  along  the  lengA  of  the  specimen  at  two  locations.  The  first, 
placed  at  the  mid-plane  free  edges  (30/-30  interface),  while  the  second,  symmetrically 
placed  at  60/60  and  60/-60  interfaces.  The  damage  was  symmetrically  stacl^  in  the  second 
configuration,  in  orcter  to  avoid  warping  due  to  curing  stresses. 

The  specimens  were  tested  under  uniaxial  loading  and  the  associated  end-twist  was 
recorded.  Hie  results  showed  a  reduction  of  18.5%  from  the  undamaged  state  for  mid-  ^ 
plane  delamination  and  areduction  of  13.1  %  for  off  mid-plane  delamination  [7i]. 

Analysis  and  Failure  Prediction  in  Composite  Stiffeners: 

An  illustration  of  the  ability  of  the  developed  analytical  models  to  predict  damage 
onset  and  growth  at  the  component  level  as  well,  is  provided  in  Refs.  72-/5.  These 
include  composite  plates  and  stiffeners  subjected  to  compressive  loads.  Of  significance,  is 
the  result  of  Ref.  15  where  the  failure  sites  observed  in  tests  of  composite  stiffeners  are 
explained  based  on  the  interlaminar  stress  field  associated  with  the  postbuckled 
configuration.  This  correlation  was  shown  in  Ref.  9  for  damage  onset  inediction  in 
unsymmetrical  laminates.  A  similarity  between  failure  prediction  based  on  interlaminar 
stresses  and  strain  energy  release  rate  was  also  established 

Lessons  Learned  and  Future  Investigation: 

Among  the  findings  of  this  research  woiic,  three  lessons  have  significant  implications. 

The  first,  is  the  similarity  in  analytical  modeling  between  syrnmetrical  laminates  subjected 
to  applied  combined  loads  and  unsymmetrical  laminates  subjected  to  separate  loading.  In 
the  latter  the  combined  loading  effect  is  induced  by  the  intrinsic  coupling  associated  with 
the  unsymmetry.  The  second,  is  the  similarity  in  ^ure  predictions  between  stress  based 
models  and  strain  energy  release  rate  approaches.  Finally,  fiee-edge  delamination  can  have 
a  significant  effect  on  the  coupling  stiffness  of  a  laminate  while  its  influence  on  the  axial 
stiffoess  is  ne^dble.  The  variation  in  the  coi^ling  stiffness  depends  on  the  delamination 
site  and  can  lx  dctKted  at  a  low  level  of  loading.  This  finding  points  to  a  basic  inquiry, 
namely,  how  a  given  damage  mode  affects  the  various  anisotropic  stiffness  coefficients  of 
laminated  composites.  Also,  if  the  variation  in  a  stiffness  coefficient  could  be  directly 
correlated  to  a  given  damage  mode  and  site.  These  inquiries  provide  the  motivation  for  the 
proposed  stractures  and  materials  research  task  mi  dainage  characterization  and  analysis  of 
unsymmetrical  laixtinates. 
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External  Interactions 

A  number  of  briefings  has  been  made  to  the  Army  Aerostmctures  Directorate  at 
NASA  Langley  and  Sikorsky  Aircraft  The  sublacninate  analysis  of  laminates  subjec^ia 
bending  and  combined  loa^  was  provided  to  Dr.  Habib  RaLof  Bell  Helicower.  The 
Siil^sis~bf~ unsymmetrical  laminates  was  given  to.  Dr._RoDaldJZabora  oC  Boeing-  *■ 

CommerdiJ  Aii^Iw^lh*.'  Rodenck  Martin  of  NASA  Langley  and  Dr.  Steven  Hooper  at _ 

Wichita  State  University.  The  analysis  and  design  methodology  of  lamirates  with  optimum 
extension-twist  coupling  were  provided  to  I>r.  Ran  Kim  of  the  University  of  Dayton 
Research  Institute. 
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IV.  Flight  Mechanics  and  Controls 


Task  1.  Modern  and  Active  Control  Research  for  Rotorcraft 
Applications 


AJ.  Calise,  D.P.  Schrage  and  J.V.R.  Prasad 


Problem  Studied 

This  task  is  devoted  to  exploring  robustness  and  controller  order  reduction  issues 
related  to  the  design  of  rotorcraft  flight  control  systems,  and  active  control  of  bla^  modal 
responses  to  disturbances  and  pilot  inputs.  Tbe  bbjtvTives  are  to  develop  a  methodology  for 
designing  fixed  order  dynamic  com^nsators  in  an  optimal  output  fe^back  setting.  The 
design  approach  should  allow  for  a  dW:t  means  of  tra^g  off  performance  and  robusmess 
measures.  Applications  to  rotorcraft  flight  control,  and  active  control  oi  rotor  dynamics  and 
blade  modal  responses  have  been  investigated.  Robusmess  to  both  unstructured  uncertainty 
and  to  structured  parameter  uncertainty  have  also  been  addressrxl.  For  simultaneous  control 
of  rigid  body  and  aeroelastic  modes,  a  two  time  scale  design  approach  was  develoi^ 
based  on  Singular  Perturbatitxi  Theory.  Ntmlinear  transformation  techniques  are  also  being 
investigated  which  provide  a  means  of  develc^ing  invariant  controllers  that  give  a  desired 
response  in  all  flight  nxxles. 

Summary  olRcsBlts 

During  this  pmod  we  have  continued  our  work  to  develop  necessary  and  sufficient 
conditions  for  designing  fixed  order  dynamic  compensators  that  satisfy  an  Hee  bound,  and 
which  in  the  limit  as  a  parameter  is  reduced  will  peput  the  calculation  of  the  optimal  fixed 
order  Hoo  controller.  We  are  very  close  to  corrqrleting  the  sufficiency  part  of  the  proof,  and 
have  gained  considerable  insight  into  the  conditions  that  must  be  imposed  before 
sufficiency  can  be  denoons^rated.  Several  fairly  large  order  fli^t  control  examples  are  also 
being  developed  which  will  be  used  to  illustrate  the  metho&logy  and  benefits  of  fixed 
order  compensation.  A  paper  is  cunendy  being  prepared  for  submission  to  the  Decision 
and  Control  Conference  in  which  we  hope  to  summarize  the  developments. 

We  have  also  initiatea  a  new  effort  to  improve  the  numerical  method  used  to 
compute  an  t^timal  outyut  feedback  gain.  This  area  remains  a  major  stumbling  block  to  the 
development  of  reliable  software  tor  opdmizuig  the  gains  in  a  fixed  conqwnsator  structure. 
Two  approphes  are  currently  under  investigation.  The  first  uses  a  novel  homotopy 
method  which  starts  with  the  solution  for  full  state  feedback,  and  progress  to  the  outyut 
feedb^k  solution  by  following  the  solution  path  for  a  set  of  algebraic  equations  as  a  sc^ar 
pararx»ter  is  varied  from  zero  to  (xie.  When  the  parariKter  is  zero,  soluti(»  of  the  algebraic 
equations  is  equivalent  to  die  full  state  feedback  solution.  When  the  parameter  is  one,  the 
solution  solves  the  ouqiut  feedback  inoblem.  The  potential  advantages  to  diis  approach  are: 
(1)  it  eliminates  the  nekl  for  specifying  an  initially  stabilizing  output  feedback  gain  (needed 
as  a  starting  point  for  tbs  current  algorithm),  and  (2)  if  the  solution  path  (udiich  is  followed 
by  planning  a  quadrature)  does  not  contain  a  singularity  the  t^timal  ou^ut  feedback 
solution  is  obtained  without  the  need  for  iteratitm.  The  second  approach  entails  modifying 
fee  current  iterative  sequential  algorithm  so  that  it  searches  txdy  over  fee  set  of  intnnally 
stabilizing  output  feedmick  gains.  A  noajor  deficiency  in  die  existing  iqiproach  is  feat  it 
often  takes  a  search  direction  which  may  bad  to  destabilizing  gain  for  a  very  small  step. 
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When  this  occurs,  the  algorithm  fails  to  find  a  minimizing  solution,  even  though  we  have 
developed  a  proof  that  for  a  sufficiently  small  step,  the  p^oimance  index  can  always  be 
reduced  The  difficulty  is  that  it  is  possible  that  die  improving  step  may  become  arbitrarily 
snoall,  and  the  resulting  improvement  may  not  be  computable  within  machine  precision.  It 
is  hoped  that  by  restricting  the  search  to  internally  stabilising  gains,  this  difficulty  will  be 
avoided.  This  restriction  entails  enforcing  a  linear  constraint  on  the  controller  gains  in  the 
search  process.  We  have  developed  a  proof  of  convergence  for  a  sequential  algorithm 
which  shows  that  given  an  initially  stabilizing  gain  which  doesn't  necessarily  satisfy  the 
constraint,  the  constraint  will  be  satisfied  afin'  a  finite  number  of  steps  and  remain  satisfied 
thereafter  in  each  iteration  of  the  search  process. 

During  this  period  we  have  considered  approximate  full  model  inversion  for 
nonlinear  controller  synthesis  using  the  feedback  linearization  technique.  The  control  terms 
that  appear  in  the  nonlinear  equations  were  approximated  in  terms  of  a  control  sensitivity 
matrix  in  order  to  simplify  the  inversion  pro^ss.  Dynamic  inversion  essentially  involves 
inversion  of  control  sensitivities  in  order  to  calculate  the  required  control  movements.  In 
the  helicopter  case,  certain  control  sensitivities  are  very  small,  e.g.,  cyclic  control 
sensitivities  on  body  axis  accelerations  are  generally  very  small.  It  is  seen,  using  the  full 
model  inversion  results,  that  the  effect  of  sn^  control  sensitivities  is  to  cause  significant 
control  chattering  and  as  the  bandwidth  of  the  controller  is  increased,  not  only  does  the 
chattering  increases  in  frequency  but  also  the  closed  loop  system  eventually  goes  unstable. 
It  is  felt  that  an  alternate  way  to  use  an  approximate  mo^l  that  includes  terms  contaiiung 
products  of  system  states  and  controls.  We  are  currently  investigating  the  alternate 
approach.  Also,  during  this  period,  we  have  conceptualized  a  general  procedure  for 
synthesizing  robust  nonlinear  controllers  for  helicoptos.  The  procedure  involves  the  use  of 
system  identification  techniques  for  obtaining  an  approximate  nonlinear  nxxlel  of  the 
helicopter  in  polynomial  form  of  prespecified  degree.  Then  api^ximate  linearization 
techniques  are  us^  to  transform  the  approximate  nonlinear  model  into  linear  domain.  The 
issue  of  describing  maximum  bounds  on  uncertainties,  diat  are  needed  for  robust  controller 
design  in  the  linear  domain,  can  be  addressed  using  fit  errors  between  flight  test  data  used 
in  the  identification  step  and  response  data  con:^>ut^  using  the  approximate  model.  Hnally, 
linear  controller  synthesis  is  performed  on  the  transformed  system  in  order  to  find  the 
feedback  laws  necessary  to  obtain  performance  in  the  presence  of  uncertainty.  A  paper 
describing  this  approach  has  been  accepted  for  presentation  at  the  1993  Annual  Ftirtm  of 
the  American  Helicopter  Society. 
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